Wiggles in the cosmic microwave background radiation: echoes from non-singular 

cyclic-inflation 
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In this paper we consider a unique model of inflation wfiere tfie universe undergoes rapid asym- 
metric oscillations, each cycle lasting ~ 10* Planck times. Over many-many cycles the space-time 
expands to mimic the standard inflationary scenario. Moreover, these rapid oscillations leave a 
distinctive periodic signature in In k in the primordial power spectrum, where k denotes the comov- 
ing scale. Although the cyclic-inflation model contains additional parameters as compared to the 
standard power-law spectrum, the improvement to the fit of the 7-year WMAP data is significant. 



I. INTRODUCTION 

Primordial inflation has been very successful in ex- 
plaining the perturbations in the cosmic microwave back- 
ground (CMB) radiation and the large scale structures 
of the universe [l|, for a recent review see 0- However 
inflation has some outstanding problems, in particular 
it doesn't encode a non-singular geodesically complete 
evolution We will present a unique singularity free 
geodesicahy complete realization of inflation in the con- 
text of cyclic cosmologies. 

In cyclic cosmological models, rather than having a sin- 
gular beginning of time, our universe can be made eternal 
in both past and future P-[l3|. However, in most cyclic 
scenarios the effort has been to produce primodial fluctu- 
ations within a single long cycle ^. Although there have 
been progress jl5l4l7l | , this has proved rather challenging 
in comparison to the success which inflation enjoys in ex- 
plaining the observed near scale invariant perturbations 
in the CMB. 

In this paper we provide a simple alternative to the 
standard inflation and cyclic universe scenarios in the 
form of a cyclic-inflation model which tries to incorpo- 
rate the successes of both; our model includes a non- 
singular cyclic phase of evolution where in every cycle 
the universe contracts a little less than it expands lead- 
ing to an overall growth. In fact, over many-many cycles 
the space-time resembles that of inflation. Thus, in close 
analogy with inflation, we can explain how the seed per- 
turbations generated at much higher energy densities can 
be stretched to the observable scales, and why the spec- 
trum is nearly scale free. Additionally, the model leaves 
distinct signatures of the rapid oscillations the universe 
undergoes by modifying the power spectrum with peri- 
odic signatures. Last but not the least, it turns out that 
the cyclic-inflationary phase requires a negative potential 



energy, but the universe can gracefully exit to a positive 
potential region marking the end of the inflationary phase 
and the onset of standard radiation dominated era. Thus 
our model may provide a way of reaching a positive en- 
ergy vacuum from a plethora of negative energy vacua 
in the string landscape (l8| . Finally, Tolman's entropy 
problem (which is equivalent to the problem of geodesic 
incompleteness in our model) can be naturally addressed 
by including a pre-inflationary emergent phase where the 
scale-factor starts to oscillate periodically as we approach 
past infinity, the size of the universe never becomes van- 
ishingly small [loj . 

Let us consider a simple cyclic inflation model, where 
the universe is mostly dominated by radiation, and the 
cycles are (approximately) periodic in energy densities. 
This follows if firstly, we assume that quantum gravi- 
tational effects trigger a bounce whenever some critical 
Planckian energy density is reached. Secondly, we need a 
— ve cosmological constant (CC), —A, which ensures that 
the universe turns around and starts to re-collapse once 
the radiation energy density dilutes and becomes equal 
to A. Thus contrary to common expectations, in the 
presence of matter the universe does not get stuck in an 
AdS vacua but rather starts to cycle. These cy- 

cles are typically short, the time period, r — oMp/^/X, 
where a - 0(1) and Mp = 2.4 x 10^* GeV [l^- We 
shall show that in order to obtain the correct amplitude 
of CMB fluctuations wc will require A to be close to the 
conventional string/GUT scale, A^/^ ^ 10 



^Mp, so that 



^ Recently, there have been attempts to calculate how perturba- 
tions can evolve through various cycles |14| . 



Now, provided there is exchange of energy between ra- 
diation and some other forms of matter, then, as a nat- 
ural consequence of the second law of thermodynamics, 
one expects the cycles to be asymmetric. The total en- 
tropy in the universe increases monotonically, and by the 
same factor in every cycle: Sn+i/Sn = 1 + 3k. Since 
entropy is proportional to the volume this means that if, 
for instance, we track the scale factor at the bounce point 
of consecutive cycles then 

a„+i/a„ = 1 + K for K <C 1 (1) 

While the above scenario can be realized in many differ- 
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ent ways, here we will consider a simple toy model with 
two species, massless radiation and some massive particle 
which interact with each other. It is clear from ([T]) that, 
over many asymmetric cycles the evolution of the uni- 
verse looks very similar to that of ordinary inflation with 
an average Hubble expansion rate H^m ~ k/t. We will 
see later that this "cyclic inflationary" phase can indeed 
address the usual cosmological puzzles such as isotropy, 
horizon, flatness and homogeneity. 



Last Bounce 




Inflation 



Exit 



FIG. 1: Typical potentials: In the positive energy side three 
different possibilities are consistent with our model. 

What about the spectrum of the primordial fluctua- 
tions? To match the COBE normalization, the power- 
spectrum associated with metric fluctuations must be 
given by: 7-"$ ~ 10~^°. Now in general, since mat- 
ter couples very weekly to gravity, in the sub-Hubblc 
phase (when the wavelength of a given comoving mode 
is smaller than the cosmological time-scale), when the 
metric fluctuations are generated from the matter fluc- 
tuations the amplitude is suppressed by the Planck scale. 
Typically we have 

r-i, cx k^^l - p/M^ - 10-1° . (2) 

Once the wavelength becomes larger than the cosmolog- 
ical time scale, the metric fluctuations effectively freeze 
at the value of p when the particular mode crosses the 
Hubble radius. This intuitive picture will essentially let 
us argue why the perturbations in our model will have 
a near scale invariant spectrum with a distinctive peri- 
odic feature (in In A:) that, in fact, provides a significant 
improvement in the WMAP 7-yr fit, and hence may be 
detectable in the future experiments. The two most im- 
portant parameters governing the physics are A and k. 
While the former determines the amplitude of fiuctua- 
tions, the latter characterizes the wiggles on top of the 
near scale-invariant spectrum. 

Finally, let us discuss the graceful exit problem in this 
infiationary scenario. If we are stuck in a — ve CC, then 
the above infiationary phase persists forever and one can 
never obtain an universe like ours. Fortunately, one can 
exit the inflationary phase if instead of a negative cosmo- 
logical constant we have a dynamical scalar field whose 
potential interpolates between a negative and a posi- 
tive cosmological constant as depicted in Fig. [T] Since 
V{(j)) — ^ —A as (/) — >■ oo, we can realize the inflationary 
phase, but the scalar held keeps rolling towards smaller 



(f) and eventually there comes a (last) cycle, in this last 
contraction phase the scalar fleld gains enough energy 
to zoom through the minimum and reach the -|-ve CC 
phase. 

The paper is organized as follows: In the following 
section [III we describe our toy model and the back- 
ground evolution which mimics the inflationary space- 
time. Next, in section iHll we explain why we expect to 
get a nearly scale-invariant spectrum with characteristic 
small wiggles in our model. We also try to flt our model 
with the WMAP data and report our findings. In II VI we 
end with a discussion of the standard cosmological puz- 
zles in the context of our model and an outlook towards 
future research directions. 



II. THE MODEL AND BACKGROUND 
EVOLUTION 

A. The Cycles 

Let us start by looking at a simple toy model where 
we have a negative cosmological constant, and the "mat- 
ter content" of the universe consists of a single non- 
relativistic species along with relativistic degrees of free- 
dom. The scenario we will present here is very similar to 
what was considered in [l^, except that we will have a 
radiation dominated universe, where as in [T2j the domi- 
nant contribution to energy density always came from the 
non-relativistic species. The advantage of having a radia- 
tion dominated universe is that the model can then easily 
avoid the Black-hole over-production problem, common 
in cyclic models, as the Jeans length is always large and 
comparable to the Hubble length. Most of the analysis in 
the present radiation dominated scenario closely parallels 
what was carried out in [13] , only some numerical factors 
change. Hence, we will keep this part of the discussion 
brief, and refer the readers to the appendix for details. 

To obtain the cyclic evolotuion, we are going to make a 
couple of assumptions: Firstly, we will assume that dur- 
ing contraction if a critical Planckian energy density, p^, 
is reached, the universe bounces back nonsingularly to 
a phase of expansion ^. We will see later that the de- 
tails of the mechanism of the bounce is not particularly 
important for our model. Secondly we are going to de- 
mand that the relativistic and the non-relativistic species 
remain in thermal equilibrium above a certain critical 
temperature, Tc, but below this temperature the massive 
non-relativistic degrees of freedom fall out of equilibrium, 
and consequently when they decay into radiation thermal 
entropy is generated. We will see shortly that this is the 



^ Although several different bouncing models have been considered 
in the literature, in most models a bounce occurs when the energy 
density reaches close to the Planck density, see for instance [a, 
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crucial ingredient that makes our cycles slightly asym- 
metric and causes the universe to "effectively" inflate. 
In the non-thermal phase the Hubble equation reads: 



rp4 



3A/2 \ a' 



A 

T4 



(3) 



where —A = —M^ is the negative cosmological constant, 
and the fi's are related to the energy densities via 



rp4 



and pr = 



4 



(4) 



For definiteness, we consider a compact universe ^ with 
a volume V = T~^a^. Let us define the ratio of the equi- 
librium energy densities of the non relativistic (massive) 
and relativistic (massless) species at the transition point, 
T = Tc, to be given by 



(5) 



We will here be interested in scenarios where /i ^ 1, so 
that to a good approximation the universe is always radi- 
ation dominated. In order to understand the dynamics, 
it is instructive to first look into the evolution when non- 
relativistic and relativistic species are non-interacting. 
Then ft's would just be constants yielding a Friedmann- 
Robertson- Walker (FRW) cosmology, and we will have a 
cyclic universe scenario. Let us start tracking the evo- 
lution just as the universe enters the non-thermal phase 
at the transition temperature, Tc, during the expanding 
branch. As the universe expands further, the tempera- 
ture decreases until T ^ Ai at which point the universe 
turns around due to the presence of the negative cosmo- 
logical constant. As it contracts the temperature rises. 
Once the temperature crosses Tj,, the universe enters the 
thermal (still contracting) phase. Once Planckian den- 
sities are reached, according to our prior assumption, a 
quantum bounce transitions the universe back to the ex- 
panding branch. Once more the temperature dilutes and 
as it falls below Tc, the next cycle begins. One can com- 
pute the time period of these cycles just by integrating 
the Hubble equation ([3]). In appendix |V] we have ob- 
tained the expression in the approximation where wc ne- 
glect non-relativistic matter contribution to the Hubble 
equation ^ (see also [l^l) ^: 



2VA 



(6) 



Importantly, this is a constant and does not change 
from cycle to cycle as it does not depend on fi^ , (As 



^ For an open or a flat universe one just has to rephrase all the 
arguments in terms of entropy density rather than the total en- 
tropy and volume of the universe. 

We have ignored the time spent in the thermal bounce phase, 
because approximately it is given by: T(, ~ Mp/T^, which is 
much shorter than t as long as T* 3> A. 



we will see shortly, the fi's will increase from cycle to 
cycle with the gradual expansion of the universe.). An- 
other way to think about this is that although the cycles 
are asymmetric in scale factor, it is periodic in energy 
densities, and it is the various energy densities that the 
Hubble rate is sensitive to. 



B. Energy Exchange & Inflation 

Let us now turn our attention to the overall growth of 
the scale-factor, which is best expressed in terms of en- 
tropy growth from cycle to cycle. The amount by which 
the entropy grows as matter gets converted into radiation 
depends on the details of the microphysical processes in- 
volved, but for us the only thing that is important to 
realize is that since the different p's and H(t) are peri- 
odic functions of time, this growth in entropy density will 
not change from cycle to cycle. For the purpose of illus- 
tration in appendix |Vl we considered a very simple sce- 
nario where in the non-thermal phase the non-relativistic 
species can simply decay into the radiative degrees. If the 
decay time is much larger than the time period of the cy- 
cle. Ft <^ 1, one can analytically compute the increase 
in the scale-factor in a given cycle: 



1.6pTMpTcg 

A4 



1/4 



(7) 



As expected the increase in entropy depends on the var- 
ious mass scales involved, the parameter n which con- 
trols how much of non-relativistic matter is present in 
the fluid, and F which determines how fast the massive 
particles decay into radiation. It is easy to see that by 
choosing appropriate parameter values one can make this 
number small. 

Now, in a more realistic paradigm one will have several 
species involved. A stringy model could involve a ther- 
mal Hagedorn phase near the bounce when the massive 
string states are kep t in thermal equilibrium with the 
massless degrees [l0|, [il, but below some critical tem- 
perature when the relevant scattering cross-sections be- 
come too small, they would fall out of equilibrium. If 
these modes consequently decay into radiation, entropy 
would be produced in a manner very similar to what 
we discussed above. After the turnaround, as the uni- 
verse starts to contract and the energy density increases, 
the massive states can be replenished from radiation by 
scattering processes. In this case the parameter k would 
capture stringy physics involving thermodynamics, scat- 
tering cross-sections and decay rates. 

Finally, let us comment briefly on the graceful exit 
mechanism. Once we replace the negative CC with a 
negative potential of the form depicted in figure [1] the 
universe can exit the cyclic inflationary mechanism. This 
is possible because, if the slope of the potential suddenly 
becomes steep the total scalar energy density can be- 
come positive during a single contraction phase. Once 
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the total scalar energy density is positive, since the en- 
ergy density can only increase in the contracting phase, 
the scalar field cannot turn back in the negative poten- 
tial region (at the turning point this would imply nega- 
tive total energy) . This is a well known result, see for 
instance [l2l . [l9j . The universe bounces one last time 
when the energy density reaches ^ pb- As long as the 
bounce occurs when the scalar field is already in, or is 
"sufficiently near" , the positive part of the potential, the 
present universe will emerge dynamically with a posi- 
tive cosmological constant after the graceful exit from 
the cycling inflationary phase. Note, that the universe 
cannot turnaround any more as the scalar energy den- 
sity is no longer negative. Moreover, after the bounce 
since the kinetic energy of the scalar field starts to red- 
shift as a~^, even if it dominates matter/radiation, it 
will quickly become subdominant ensuring the entry into 
a matter/radiation dominated universe. 

Approximate calculations corroborating the above ar- 
gument was presented in p^ . and we are currently per- 
forming a detailed numerical exploration (23j of the entire 
parameter space to determine whether any fine-tuning is 
required for the mechanism to succeed. 



metric power spectrum is scale- free: $fc — > k^'^l"^ ^ . Ac- 
cording to our previous arguments therefore, as pertur- 
bations become super-Hubble (t — )■ oo and A — >■ oo) this 
scale- invariance should be preserved, modulo periodic de- 
viations: 

= ^ T'i.anfc -f ln(l + K)) = T'^anfc) 

(9) 

Since K is typically a small number in our model, we 
expect to have a near scale-invariant spectrum. 




fl-l ■ 1 ■ 1 ■ 1 ■ r 

12 3 4 



III. CMB SPECTRUM 

A. Scale-invariant envelope with Periodic 
Modulations 

Let us first provide a general argument as to why 
we expect to obtain a nearly scale-invariant spectrum 
in our model. We are interested in tracking As 
the scale-factor evolves, so does the physical wavelength, 
Afc(t). Now the evolution of $fe at a given physical wave- 
length, A, depends on the background environment that 
the mode "sees" , i. e. it depends on the different energy 
densities involved, {pi{ty\, where i labels the various fiuid 
components. Thus another way of looking at the prob- 
lem is to realize that the evolution of $fc depends on the 
curves {yOi(A)}, and these curves will in general be differ- 
ent for different comoving modes. For simplicity, let us 
first consider the case when we a — ve CC (then cycles 
are precisely periodic in energy densities), and look at 
two modes which are related by k, the factor by which 
the universe expands in a given cycle: 

fc' = (1 + K)k ^ (1 + K)Afc' = Afe (8) 

The curves {pi{X)} for these two modes are completely 
identical, the fc'-mode lagging behind a cycle as compared 
to the fc-mode. Therefore as one evolves from A 0, the 
sub-Hubble phase, to A — >■ oo, the super-Hubble phase 
where the fluctuations freeze, these two modes should 
change by the same factor. If we use the traditional 
Bunch-Davis vacuum initial conditions at < — >■ — oo (or 
A 0) like in standard inflation, then the sub-Hubble 



FIG. 2: The Last Exit: The curve in the red corresponds to 
the cosmological time-scale, while the green and the yellow 
curves represent the physical wavelengths of two different co- 
moving modes. The two modes are initially in a Mixed phase 
but then they make their last exits after which they evolve 
as super Hubble modes. While the green wavelength exits in 
the third cycle, the yellow curve has to wait for another cycle 
to make it's last exit. The two modes experiences identical 
background cosmology, the yellow curve only lagging behind 
by a cycle as compared to the green curve. 

Let us actually try to calculate the spectrum by mak- 
ing a few simple assumptions. Firstly, we are going to 
assume that as long as wavelengths of fluctuations are 
smaller than the cosmological time scale their power spec- 
trum is determined by the energy density of the fluid, vis 
a vie ^ . Secondly, we will assume that once the modes 
become super-Hubble they stop evolving and freeze out. 
These are indeed what one flnds in GR and as we will see, 
in our scenario the freeze-out occurs near the turnaround, 
away from the non-GR bounce. Hence we expect these 
assumptions to hold true but we do want to caution the 
readers about possible modifications coming from the 
non-GR bouncing phase. 

Having said that, let us go ahead and try to calculate 
the amplitude of perturbations, which essentially boils 



Bunch-Davis vacuum quantum fluctuations are implemented 
on the canonical Mukhanov variable Vf^ which is nevertheless 
straight-forwardly related to In particular the quantum vac- 
uum fluctuations, Vf^ ~ k~^f'^ yields ~ k~^/^, see for a re- 
view [2||. 
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down to calculating the amplitude at the time of the 
"last exit" . To determine when this last exit occurs let's 
first look at the cosmological time scale. In the usual 
monotonically expanding universe this is just given by 
the Hubble radius. However, during the bounce the time 
scale, Tb is determined by the bounce energy density, pi,. 
Typically 



n 



(10) 



where /3 ~ C'(l) constant. The cosmological time 
scale thus approximately behaves as follows: During the 
bounce phase it stays a constant and is given by t;,, 
(fTO|) . After the bounce, radiation dominated GR takes 
over, the cosmological scale just corresponds to the Hub- 
ble scale, and starts to increase. Near the turnaround 
phase, it reaches it's maximum value which is essentially 
the time period of the cycle determined by cosmological 
constant, see (jG)). After the turnaround, the cosmologi- 
cal scale decreases back to it's minimum value near the 
bounce, and this cycle repeats itself, see fig. [21 

In the mean time, the physical scale corresponding a 
particular comoving perturbation keeps oscillating with 
the cycles of the universe, but because of the asymme- 
try gradually expands. Thus, there is an initial "pure 
sub-Hubble" phase where the fluctuations are always 
smaller than the cosmological scale. Then there appears 
a "mixed phase" where the fluctuations go in and out 
of the "Hubble radius" . Finally, there comes a last cy- 
cle where the wavelength goes out of the Hubble radius 
never to return, see figure [2l We call this the last exit. 
According to the two assumptions we have made, all we 
have to do is to calculate the amplitude of fluctuations 
of the modes at their last exit point. 

For this purpose let us focus on a single cycle. Con- 
sider modes which are super- Hubble at the bounce, i.e. 
Afe ^ Tf). After the bounce in the radiation dominated ex- 
pansion phase, the Hubble radius grows faster than the 
scale factor, and hence these modes start to re-enter the 
Hubble radius. This process continues till we reach the 
turnaround point when the Hubble radius stops grow- 
ing and therefore the modes stop re-entering. After 
the turnaround in the radiation dominated contraction 
phase, precisely the opposite happens, i.e. the modes 
start to exit the Hubble radius. Now, and this is the 
key point, since in our model the universe contracts a lit- 
tle less than it expands, those modes which are the first 
to exit the Hubble radius in the contraction phase, never 
come back inside the Hubble radius in the following cy- 
cle. They have just made their last exit after which their 
amplitudes freeze. This is illustrated in figure [21 

We are now ready to calculate the spectrum. To this 
end we need to compute the energy density of radiation 
at the time the modes exit during contraction. The ra- 
diation density is, as usual, given by 



Pb 



at 



To obtain the spectral dependence we need to use the 
Hubble crossing condition 



a/k = H' 



k cx {-ty 



(12) 



since during the radiation dominate contraction, H oc 
1 /t, and a oc Thus we have 



Pr <x (— t) ^ cx A:^ 



(13) 



If ki corresponds to the comoving wave which is the first 
one to exit the Hubble radius during contraction, then 
all the modes with values up to A:i(l -I- k) also exits the 
Mixed phase in the same cycle. Thus the power spectrum 
for the modes that exit the Mixed phase in a given cycle 
reads 



A 
\k 



for fci < fc < (1 + K)ki (14) 



(11) 



7 being an 0(1) parameter relating the Power-spectrum 
to the energy density of the fluid; we note that since 
the modes exit near the turnaround the energy density 
in radiation which controls the amplitude ^ must be 
close to A. Now, once the wave number increases by a 
factor (1 -|- k), the fluctuations can no longer exit the 
Mixed phase, because they re-enter the Hubble radius 
during the following expansion. They therefore, must 
wait for the next cycle for their last exit, see flgure [H 
and consequently the power spectrum repeats. 

A few remarks are in order. Although our arguments 
are quite general and certainly valid in GR, it is pos- 
sible that the Quantum bounce regime may modify the 
analysis. To perform a complete treatment, one would 
need to work with a concrete model for the bounce, and 
then track the perturbations (most likely numerically) 
as they evolve through various cycles. We are currently 
engaged in such an analysis [23j . It is worth point out 
though that, given that the short wavelengths tend to de- 
couple from cosmology it seems a reasonable assumption 
that they won't feel the effect of the bounce and hence 
their fluctuations will continue to be given by the stan- 
dard expression [21 Also, it is hard to imagine that the 
super-Hubble constant mode in GR will disappear in the 
more general theories. See for instance [l^l, for a specific 
higher derivative example where it was found that the 
super-Hubble modes indeed freeze out. 



B. Spectral Tilt 

The fact that the energy density of the scalar field dur- 
ing the inflationary phase is not a constant but slowly 
evolves implies that the spectral envelope (without the 
periodic wiggles) is not going to be precisely scale- 
invariant but have a tilt. Now we are working on the 
assumption that the amplitude of oscillation freezes when 
the physical wavelength of a given comoving mode equals 
the cosmological time scale. This in turn depends at 



what scale the turnaround occurs. In our scenario, the 
turnaround occurs when the positive radiation density 
cancels the negative scalar field density. Thus, 7-"$ ex Ip^l, 
and the power spectrum will change slowly as </> evolves. 

Now, if V'{(l>) = 0, 4> would just get Hubble damped 
and freeze. However, for a non-zero V'{(t>) it will evolve 
in a manner very similar to the inflationary slow roll: 
3Hav(}) = -V'{(j)), so that 



(15) 



i8y(0)K2 

Depending upon the parameters of the potential it is pos- 
sible for p{(t)) to cither increase or decrease as (j) evolves 
resulting in a red or a blue tilt respectively. For in- 
stance, if our potential in the inflating region is modeled 

as = -A(l + e-'^'^/^^p), with e~'"t'^'^'^p < 1, then 

one finds 

P0 « -A[l + e-"l'/^^^ - ((^2^718^2)6-2'^'^/^^^] (16) 

Since k is a small parameter one can check that it is 
possible for the third term to dominate over the second 
so that we obtain a red tilt. We plan to investigate the 
detail dynamics in future, the important point for now 
is that, in general we will have a non-zero spectral tilt 
which can be either red or blue. Thus we expect the 
power spectrum to have the following approximate form: 

-1 



Pb 



f{k) 



(17) 



where f{k) is defined as 
fik) 

& f{k{l + K)) 



k 

fik) 



for ki < k < {1 + K)ki 



(18) 



In particular, a distinctive feature of these modulations 
is that / varies from 1 to 1 4k as the wave number goes 
from fci to (1 4- K)ki. In other words, the amplitude of 
fluctuations. A/ = 2k, is related to the period of fluctua- 
tions, A In A: = K in a simple way. This is important since 
previously very different physical models have been con- 
sidered which give similar oscillatory modulations of the 
spectrum, but the amplitude and the period arc typically 
two independent parameters [23 |. 



C. Fitting WMAP-7yr data 

To estimate the effect of the wiggles as a good approx- 
imation to the discontinuous wiggles (jlSp . one can con- 
sider a smooth oscillatory modulation since it has almost 
the same shape and oscillatory behavior and is also more 
appropriate to be used to calculate the angular power 
spectrum using CAMB or COSMOMC. We can there- 
fore use a power spectrum of the form: 



A 



1 + 2k cos 



2tt , k 




(19) 



FIG. 3: Resultant best fit (top panel) and (bottom 
panel) from the cyclic inflation in comparision with 7-year 
WMAP binned data. 



where fco is an arbitrary normalization scale. The main 
additional parameter as compared to the standard pri- 
mordial spectrum is k < 1/2, although we also have an 
additional phase factor in 6' G [0, 27t]. 

In order to see how well the cyclic inflation flts the 
CMB data, we have done some coarse sampling in the pa- 
rameter space. For baryonic abundance; ^li,h^ — 0.0228, 
cold dark matter abundance, ildmh^ = 0.1156, h — 0.70, 
the optical depth, r = 0.082, and parameters of the 
primordial spectrum of ris = 0.963, k = 0.0455 and 
9 = 2.54:6{radian), we found a very good likelihood of 
~2ln{L) = -7470.5 to 7-year WMAP data [l|. 

In Fig. [HI we see the resultant best fit angular power 
spectrum [tt and te) from the cyclic model of infiation. 
In comparison with the best fit power-law form of the 
primordial spectrum and power-law with running spec- 
tral tilt model, our model can improve the likelihood by 
Ax^ = —3.9 and Ax^ — —2.7 respectively. It is in- 
teresting that in comparison with the phenomenological 
model of power-law with running our model has a distinc- 
tively better likelihood, although the models have similar 
number of degrees of freedom. Note here that we have 
assumed a fixed value of fco = 0.05Mpc-^ in Eq. 
In comparison with the simple power-law form of the 
primordial spectrum, considering the additional free pa- 
rameters in our theory, the improvement in the likelihood 
is not as big as to claim that our model is preferred by 
the data [25|, but at the same time, the improvement in 
the likelihood is not insignificant. Allowing tensor modes 
in the analysis will not improve the best fit likelihood and 
we expect a very small value of r (tensor to scaler ratio). 
We will leave a proper cosmological parameter estima- 
tion using CMB and other cosmological data for future 
publication, but the derived likelihood and the best fit 
parameters already suggests that this model should pro- 
vide a good fit. 
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IV. THE STANDARD COSMOLOGICAL 
PUZZLES AND FUTURE OUTLOOK 

Before we conclude, let us look at the problem of 
growth of inhomogeneity/overproduction of black holes. 
The matter fluctuations, 5 = 5p/ p can only grow as long 
as their wavelengths are larger than the Jean's length, 
Aj, given by: Aj ^ CgM^j ^ where = dp/ dp, is the 
sound velocity square. Now, in our scenario the cycles 
are short, and most of the energy content is in radia- 
tion so that the sound speed is very close to the speed of 
light, and accordingly Aj ~ ^ r. In other words, 
the sub-Hubble fluctuations don't grow. On the other 
hand, once the wavelengths become larger than the cos- 
mological time scale, t, they become super-Hubble fluc- 
tuations and evolves according to the Poisson equation: 
5k = k'^^k/ {o? p), where is the Newtonian poten- 
tial characterizing the metric perturbations. Now in the 
super-Hubble phase, $fe becomes a constant while p oscil- 
lates between a minimum and a maximum energy density. 
Thus we have: 4 < A;^$fc/(a^Pmin) ^ k'^^kM^/ [a^y/K). 
In other words 5k falls as a^^ in the super Hubble phase 
just as in ordinary inflation. Finally, cyclic cosmolo- 
gies are notoriously plagued with Mixmaster chaotic be- 
havior as one approaches the Big Crunch/Bounce, since 
anisotropics grow as ~ . However, in our model the 
same reasoning that resolves the flatness problem, also 
solves this issue. Once the cyclic-inflationary phase is 
"activated" in a small patch of the universe the chaotic 
Mixmaster behavior is avoided in subsequent cycles as 
the universe becomes more and more isotropic due to 
the overall growth of the scale factor. 

To summarize we have presented a unique way of real- 
izing inflation within cyclic universe scenario, where the 
particle trajectory is geodesically complete in either past 
and future. Although in each cycle the universe expands 
only a little bit, one can obtain a large number of total 
inflationary e-foldings, A/'tot over many many cycles. As 
long as A/'tot ^ 60 it is clear that we can explain the hori- 
zon and the flatness problems. Furthermore, every cycle 
leave their quantum imprint in the CMB with potentially 
observable wiggles determined by the enhancement of the 
scale factor during every oscillation. 

We should however emphasize that there are still sev- 
eral aspects of the model which needs a more compre- 
hensive numerical investigation. For instance, what are 
the conditions necessary to maintain thermal equilibrium 
during the bounce phase? (In fact, all we really require is 
that thermal equilibrium be established at some point of 
time in every cycle, that is sufflcient to guarantee the pe- 
riodicity in energy density.) What about the graceful exit 
problem. It is clear that whether we can transition from 
the inflationary phase to the positive potential region de- 
pends on various parameters of the model, do these pa- 
rameters need to fine-tuned? Finally, what effect docs 
the bounce play in the evolution of the perturbations. 
These are important questions and we are currently pur- 
suing them using numerical simulations (23| . Finally, it 



is also foreseeable that our scenario can potentially en- 
hance non-gaussianity through stochastic resonance dur- 
ing these oscillations, again an interesting future direc- 
tion to pursue. 

AM and AS acknowledge the support of the EU Net- 
work UniverscNet (MRTN-CT-2006-035863) and thanks 
Jan Hamann and Xingang Chen for useful discussions. 



V. APPENDIX: r AND AS' 

Here we calculate the approximate time period and the 
increase in entropy in a given cycle. We will calculate 
this under the approximation that the dust component 
can be neglected as /i <C 1 and treat as a constant as 
the amount of matter decay in a given cycle is negligible. 

To calculate the time period we start by rewriting the 
Hubble equation: 



dt 



da 



\/3AIpda 



so that the time period is given by 

_ 2\/3A/p da 



Oc a 



(20) 



(21) 



where we have neglected the duration in the thermal 
bounce phase as it is going to be much shorter than the 
above integral. Now the turnaround scale factor is given 
by 



ax 



A 



4\ 1/4 



(22) 



Thus the above integral can be rc-expresscd as 

_ 2V3Mp da _ 2V3Mp ydy 



VA 



where e = ad aT <C 1. Thus we have 

ydy _ VSttA/, 



(23) 



^/i-y^ 2VA 



(24) 



In a similar manner we can proceed to calculate the ap- 
proximate entropy increase. Phenomenologically, energy 
exchanges can be captured by generalizing conservation 
equations for the two fluids to: 



mp,n = -Tts (25) 



which now includes an energy exchange term. We can 
easily compute the net entropy increase [l^ : 



a^s 



3brTc 



4p 



1/4 



a^s 



(26) 
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where we have used the usual expressions for the thermo- 
dynamic entropies associated with matter and radiation: 



1^ - lol/4j^3/4 



M 



(27) 



Here M « Tc corresponds to the mass of the non- 
rclativistic particles, and in our convention, g — 
(7r^/30)(7*, where is the number of "effective" mass- 
less degrees of freedom. 

The growth of the entropy in a given cycle will ob- 
viously depend on the phenomenological function s. For 
the simplest case, which captures the physics of the decay 
of the non-relativistic species into radiation, s is given by 



r4 



(28) 



r being the decay rate. One can now obtain the entropy 
increase in a given cycle by simultaneously solving for 
the Hubble equation ([3]), the continuity equations ([^5)) 
and entropy growth pG]) . For some special cases this can 
also be done analytically. For the purpose of illustration 
we are going to look at the case when the decay time is 
much larger than the time period of the cycle, Fr ^ 1. 
This means that basically the Jl's change very little over a 
given cycle and for the purpose of estimating the entropy 
increase from (P7)) we can treat them to be constants. We 



can approximate e' 



-rt 



(IMl), ©, (Ell) and Q we find 



1 in the expression for S. Using 



dS_ 
da 



(t-) 



(29) 



Thus the increase in entropy in a given cycle can be cal- 
culated as 



AS*: 



2ViTn„,Mp da (ir " 1 



2^/3F^7,„AJp /"'^yHEZl 



dy- 



^/A^/c Jo \/l - y 



Now, we know that at the transition temperature the 
energy density in radiation must be given by 



1/4 



(30) 



Using the expression for ay then we find 



= — = (31) 
Using the expressions for energy density and entropy we 
can also easily relate VLm with the total entropy S: 



(32) 



Thus we finally have 



AS^'-^^^f'dy. 



2v^2/c 



l.e^iTSMpTcg^ 



(33) 
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